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Abstract. - We bring attention to the fact that in order to understand existing data on neu- 
trino oscillations, and to design future experiments, it is imperative to appreciate the role of 
quantum entanglement. Once this is accounted for, the resulting energy-momentum conserving 
phenomenology requires a single new parameter related to disentanglement of a neutrino from its 
partners. This parameter may not be CP symmetric. We illustrate the new ideas, with potentially 
measurable effects, in the context of a novel experiment recently proposed by Gavrin, Gorbachev, 
Veretenkin, and Cleveland. The strongest impact of our ideas is on the resolution of various 
anomalies in neutrino oscillations and on neutrino propagation in astrophysical environments. 
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Introduction. — The experimental confirmation of 
neutrino oscillations has not only resolved the problem of 
the solar neutrino anomaly [l], shed light on atmospheric 
neutrino data [2], inspired various reactor and accelera- 
tor based experiments, but it also has implications for 
a whole range of phenomena, such as the evolution and 
explosions of type II supernovae [3 -12 . Indeed, impli- 



cations as far reaching as the interface between general 
relativity (GR) and the quantum realm have become ac- 
cessible by the discovery that neutrino oscillations provide 
a set of flavour oscillation clocks that redshift precisely as 
required by GR 
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These insights can, at least in 
principle, be tested by precision experiments that comple- 
ment the experiments on gravitationally induced effects in 
neutron 



24 25 and atom interferometry 26 



Underlying all these successes and investigations is the 
Pontecorvo assumption that a neutrino produced in a 
weak interaction mediated process may be considered as 
a linear superposition of different mass eigenstates. 



(1) 



where the index £ G {e, /i, r} represents the three flavours, 
the index i g {1,2,3} spans the three mass eigenstates, 
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and the 3x3 unitary matrix U defines and distinguishes 
the three flavours. The CP-conjugated states are identified 
with the antincutrino, and are given by 



(2) 



Indeed such a phenomenological ansatz is able to provide 
a good broad brush understanding for much of the existing 
data. 

However, theoretically, such a highly successful ansatz 
cannot be the complete description of the kinematic struc- 
ture of neutrino oscillations; a quiet reflection reveals that 
it violates conservation of the energy-momentum four vec- 
tor. This is something that is not allowed for weak in- 
teraction mediated processes which are known to be in- 
variant under spacetime translations. This fact was flrst 
realised by Goldman in 1996 27 , and a while later by 



Nauenberg [28j. Apart from some notable exceptions, the 
physics community seems to have neglected this remark- 
able observation. A small set of other papers devoted 
to the subject of energy-momentum conservation in neu- 
trino oscillations appear to have created more heat than 
light. However, the publication of a recent paper by Co- 
hen, Glashow, and Ligeti has now made it unambiguously 
clear that the conventional understanding of neutrino os- 
cillations is not entirely tenable and that the suggested 
revision must be incorporated in any analysis of neutrino 
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oscillations 29 . This revision makes a quantum entangle- 



ment of neutrinos with the accompanying decay partners 
a necessary and unavoidable element of the physics of neu- 
trino oscillations. This has fundamental consequences for 
the correct interpretation of the relevant data. 

The focus in Ref. |29j was entirely on Raghavan's proposed 
experiment [30| and on the variations in decay times ob- 
served in the GSI experiment fSlI . Here we show that im- 
plications of the new insights extend far beyond their orig- 
inal scope. Towards this end we first review the energy- 
momentum conserving neutrino oscillation formalism in a 
concrete setting. We then argue how it enters and revises 
our understanding of neutrino oscillations and its existing 
puzzles. Once this is appreciated, a whole range of new 
questions arise. Answers to these questions may affect 
not only the design and interpretations of terrestrial ex- 
periments but it may also lead to a better understanding 
of the neutrino-induced — with quantum entanglement a 
necessary ingredient — astrophysical processes. We briefly 
discuss these aspects. 



Quantum entanglement and neutrino oscilla- 
tions. — We are inspired by the following observation. 
Whenever a mass eigenstate decays into a set of other 
mass eigenstates the relevant conservation laws may in- 
duce a quantum entanglement. For instance, in the EPR- 
like decay of a spin zero particle at rest into two spin one 
half particles (in the singlet state) , the conservation of an- 
gular momentum forces the spin projections of the decay 
products to be entangled. Such an entanglement is quite 



robust 32 and it is destroyed, e.g., when the entangled 



attribute is subjected to a measurement or an interaction. 

Neutrino identification puzzle and neutrino oscillations. 

The just mentioned EPR entanglement is induced by 
the invariance of the underlying interactions under rota- 
tions. For neutrinos, the operative interaction is the weak 
interaction of the standard model. It is invariant under 
space, as well as time, translations. Therefore, with appro- 
priate caveats and parenthetic remarks about the width of 
wave functions, etc., a neutrino produced in the decay of 
a mass eigenstate (say a pion) — which we assume is an 
eigenstate of both energy and momentum — must also be 
in an eigenstate of energy and momentum. However, since 
each of the in Eq. ([!]) has different mass, Ivi) cannot si- 
multaneously be an eigenstate of energy and momentum]^ 
This fact forbids of Eq. ([I]) to be identified with a neu- 
trino produced in a weak interaction mediated process 
cf. discussion in Sec. 2]. We shall call this circumstance 
the neutrino identification puzzle. This concern was first 



emphasised by Cohen, Glashow, and Ligeti 29 . 



To be concrete, consider the production of and 



^Same remarks apply to \i>i) and Eq. (jSjl. 



the following CP conjugated processed 
7r+ — >- + ^+ 

L 



+ i^c + 



and 



L 



(3a) 



(3b) 



The neutrino identification puzzle is the assertion that 
neutrinos and antineutrinos that are thus produced can- 
not be linear superpositions of different mass eigenstates 
(in the sense of equations ([T]) and ([2])). 

The energy-momentum four vector conserving correct 

and 
and 
- re- 



identifications, as was first argued by Goldman [27 
later emphasised by Nauenberg on the one hand 



Cohen, Glashow, and Ligeti on the other hand 29 
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quire quantum entanglement between the decay products. 
For the above-considered tt^ decay at rest, the neutrino- 
muon entanglement becomes manifest by re-writing (in 
fact, by correcting) v^^ 
3b) as 



- /i+ and i/^ - 



- /i in equations ( 3a 



and 



(4a) 



(4b) 



L/kf +mj, ki, mj) | 



Si, -ki 



The notational details are adapted from Ref. [27j |^ These 
identifications are consistent with the results of refer- 
ences 28 29 ). The energy- momentum conserving gen- 
eralisation to the case when tt* are not at rest is straight- 
forward. It can be found in Ref. |29|. 

Thus, the resolution of the neutrino identification puz- 
zle resides in the replacement of the Pontecorvo sugges- 
tion contained in equations (l][2l by the Goldman sug- 
gested source-dependent identifications such as the ones 
displayed above. It is worth noting that the original title 
of Ref. [27J was 'Source Dependence of Neutrino Oscil- 
lations'. It was only in the published version, roughly a 
decade and half later, that it changed to 'Neutrino Oscilla- 
tions and Energy- Momentum Conservation'. The energy- 
momentum conserving definition of neutrinos, and their 
quantum entanglement, is necessarily source dependent. 
For instance, the produced in a 7r+ decay has a different 
quantum entanglement than a created in the process 



or those obtained from e 



^We note the decay of /x^ explicitly for later reference below. For 
simplicity, we assume tt^ at rest. 

^Here, k; represents the momentum associated with ith mass 
eigenstate, \vi), and is the mass of the ith mass eigenstate. 
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For the discussion that follows we shall call states of neu- 
trinos and antineutrinos represented by equations (l][2) 
as Pontecorvo states while the states of the type rep- 
resented by equations ( 4apb ) will be called Goldman 
states. This allows the discussion so far to be summarised 
as: weak interactions produce Goldman states, not Pon- 
tecorvo states]^ This is a simple, and simultaneously 
profound, resolution of the neutrino identification puzzle. 
Further, Cohen, Glashow, and Ligeti (CGL) [29]) show 
that flavour oscillations do not occur for the Goldman 
states. For flavour oscillations to occur a projection to 
the Pontecorvo states must be made. They arrive at the 
following conclusion (quoted below, verbatim) 

cgl-1 When the decay products of an initial state 
of well-defined momentum evolve without 
further interaction no oscillation phenom- 
ena appear. 

cgl-2 To realize oscillations the neutrino mass 
eigenstates must be disentangled. 

Since these results are rather counterintuitive, an clement 
of repetition is in order: cgl-1 translates to the fact that 
weak interactions produce Goldman states, while cgl-2 is 
the statement that for oscillations to occur a projection to 
the Pontecorvo states must occur. 

This has the consequence that, for a Goldman state pro- 
duced at the laboratory time t = and projected to a 
Pontecorvo state at time t = t^, the standard neutrino os- 
cillation probability is modified. For the simple and rep- 
resentative 2x2 setting, the modified — ^ ly^ probability 
then reads (assuming ultra-relativistic neutrinos) 



P'(l/e ^ l^e) = H{cTd - r) 

+H{r - CTd) 



l-sin2(2(?)sin2 (l.27Am2[^ ""^"^ 



(5) 



where H{t) is the standard Heaviside step function, c = 
3 X 10^ m/s, r is the running position along the neutrino 
trajectory (measured in meters, r = is the creation 
region), the source-detector distance L is in meters, the 
mass-squared difference Am^ is in eV^, and the neutrino 
energy is in MeV^ For a given ensemble of neutrinos, 
Td may be a distribution function depending on the de- 
tails of a given set up. For example, if one of the decay 



*A parenthetic remark: once the Goldman suggestion is taken se- 
riously one must in fact also incorporate the entanglement induced 
by the conservation of angular momentum (weak interactions are 
also invariant under rotations). This forces the states given in equa- 
tions ( |4a]4b[ l to be replaced by their proper singlet state expres- 
sions which contain not only the left-transforming neutrino with the 
negative helicity, but also the left-transforming neutrino with the 
positive helicity. The latter may mimic some of the signatures of a 
sterile neutrino (similar remarks hold true for antineutrinos). This 
is a subject that must be treated in a subsequent study. In order to 
not raise too many issues at once the remainder of this paper takes 
equations | |4a|4b[ l without any additional modifications. 



products is short lived, may depend on the lifetime of 
the entangled partner. 

Implications for neutrino oscillations in labora- 
tory and astrophysical environments. — An imme- 
diate consequence that follows from expression ([s]) is that 
the KARMEN-LSND anomaly can be easily resolved if 
for these experiments is close to the KARMEN source- 
detector distance in natural units 
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Turning to the proposal by Gavrin, Gorbachev, 
Veretenkin, and Cleveland (GGVC) (34j, one finds an ex- 
perimental candidate for which the proposed modification 
to the neutrino oscillation probability would become di- 
rectly detectable, provided the disentanglement time Td is 
determined by the dimensions of the reentrant tube (and 



hence to the size of the ^^Cr source). Figure 3 of Ref. [34 
gives the ratio i?2 / Ri for the capture rates for the inner 
and outer zones of the GGVC proposal. Those calcula- 
tions assume the un-modified version of Eq. ^. Fi gure 1 
below details the ratio 



R2/R1 



(6) 



with R'2/R'i calculated using the modified oscillation prob- 
ability given in Eq. ([5]). In order to gain insight into the 
dependence of /3 on the disentanglement time, Figure 2 
shows a re-evaluation of Eq. ^ with Td increased by an 
order of magnitude; for details, the reader is referred to 
the figure captions. 

It is worth emphasising that the detectability of a quan- 
tum entanglement in this experiment presupposes that the 

the 



MiniBooNE 35 and the LSND experiments 36 



low capture rates in the Ga source experiments 



37 

38 , and 



antineutrino reactor anomaly 39 , all suggest existence of 
an eV^ range mass-squared difference and that the disen- 
tanglement time Td is roughly related to the source size (for 
'low-mass' source cas^. It would, however, be extremely 
helpful if an ab initio estimate of Td was made through 
the development of a theory of quantum disentanglement 
in nuclear media. 

Furthermore, it is not a priori obvious for many experi- 
ments that Td is CP symmetric. As an example, consider 
the neutrino source of the MiniBooNE experiment which 



is produced by the decay chain in equations ( 3a I and ( 3b ) . 
If the disentanglement time Td > t^, the muon lifetime, 
then Td is necessarily CP asymmetric. This is because the 
e"*" produced in the /i+ decay, annihilate with electrons in 
the beam stop, whereas the produced in the fi~ decay, 
do not. These asymmetric interactions of the entangled 
partners of neutrinos and antineutrinos, produced in de- 



^ Clearly, we arc guided here by calculations on the mean free 
path of neutrons and the known sizes of the critical mass in nuclear 
reactors /devices. 
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Fig. 1: In drawing this figure we have chosen Td '■— to = 
0.105 meters x c. This choice corresponds to the radial size 
of the ^^^Cr source that is used in the GGVC proposal. In the 
event that the measured /3 turns out to be extremely close to 
unity, one may turn the argument around and either measure, 
or place limits on, Td- To match the GGVC proposal, the plot 
corresponds to the case sin^ 26 = 0.3. 



cays of TT^, with the beam stop affect the projection of the 
Goldman states to Pontecorvo states via a buih-in asym- 
metry in Td- 

Interestingly, the MiniBooNE source-detector distance is 
1.8 /xs, in natural units. This is quite close to the muon 
lifetime ~ 2.20 /is. This observation when combined 
with the possibility that the CP-asymmetry of tj, for the 
MiniBooNE setting may not be same as for the LSND set- 
ting (where the source-detector distance is about 0.1 /is), 
could turn out to be an important ingredient in reconciling 
these two experiments. 



Fig. 2: A complementary plot to Figure 1 for a disentanglement 
time of Td — 10 X To . To match the GGVC proposal, the 
plot corresponds to the case sin^ 29 = 0.3. In our calculations 
we have taken care of the fact that 90% of the decays of the 
^'^Cr source produce 750 keV neutrinos, while 10% of the decays 
create 430 keV neutrinos. The advantage of presenting results 
in terms of /3, rather than in terms of R'2/R'i versus R2/R1, 
lies in the fact that in calculation for (3 we do not need detailed 
information about the detector efficiency, it cancels out. 



Conclusion. — We have provided a brief review of 
why simultaneous conservation of energy and momentum 
requires the decay partners of neutrinos to be entangled. 
No flavour oscillations occur until this entanglement is de- 
stroyed. This has consequences that have the potential 
to resolve the KARMEN-LSND conflict, change expecta- 
tions of the GGVC proposal, introduce non-intrinsic CP 
violation, and affect neutrino propagation and oscillations 
in astrophysical environments. 



It is thus apparent that quantum entanglement has the 
potential to affect the interpretation of all neutrino os- 
cillations experiments. Additionally, since neutrino oscil- 
lations play a pivotal role in type II supernovae explo- 
sions, and in the nucleosynthesis of light elements 40 



effects such as those implicit in the modified neutrino os- 
cillation phenomenology may have a significant impact for 
those studies. For pion-decay neutrinos the disentangle- 
ment time Td, interestingly, depends on the electromag- 
netic interactions, and not on the weak interactions, and 
is therefore expected to be "small". On the other hand, 
for neutrinos produced in the annihilations in neutron 
stars, Td may be quite "large" as the mean free path of 
these neutrinos in neutron-star media is of the order of a 
kilometer. In this latter instance, there are no entangled 
partners that are electromagnetically charged. This fact 
may dramatically alter the evolution of neutron stars and 
may require a recalculation of neutrino-oscillation effects 
in the explosion of type II supernovae. 
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